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1. ECUACIÓN DIFERENCIAL DE CONDUCCIÓN DEL CALOR 


Sea un paralelepípedo diferencial que forma parte de un medio en el que hay 
diferencia de temperatura de unos puntos a otros y en el que no hay cambio de fase ni 
movimientos macroscópicos entre sus partículas (advección). En estas condiciones la 
transferencia de calor, tanto en el interior del elemento como en sus límites, sólo 
puede tener lugar por conducción. 


La ecuación de la conservación de la energía aplicada al paralelepípedo es: 


dentra dE dy E d sale +U O bien dentra eN dy qe sale = U 
siendo: 
Oentra Potencia calorífica que entra en el elemento (sólo por conducción, según se 
ha justificado): 
Qentra=0x+0y+02 
om Potencia calorífica generada en el interior del elemento por unidad de 
volumen (reacción química, reacción nuclear, disipación eléctrica, ...). 
Q=Q""-dx-dy-dz 
donde q [W/m?] es la denominada generación volumétrica. 
Osalida Potencia calorífica que sale del elemento (sólo por conducción, según se ha 
justificado): 
Osalida=0x+dx FO y+dy FO z+dz 
UÚ Ritmo al que cambia la energía interna del elemento: 


: OT OT 
U=dm:c-=——= p:dx-dv:-dz-c-—— 
Ot e Ot 


donde dm es la masa del elemento diferencial de densidad p y calor 
específico c. 


Sustituyendo en la ecuación de conservación de la energía: 
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(a, + d, + q,)+ ("dx dy dz)- (4... 3 Uy +dy q dara )= pax dy de 


desarrollando en serie de Taylor y despreciando los términos de orden superior: 


A xrdx = ds An A, dx 
OX 


0q 
4 y+dy 7 q, + oy dy 


Uz+dz 7 q, + 4; dz 
Oz 


la ecuación de la energía queda como: 


0q 0q 0q hi OT 
z dx 4 2 dz +q dxdy dz = pc —- dx dy dz 
OX Oy iá Oz A iS ó Ot á 


Dado que Q,, y y q, se transmiten por conducción, en aplicación de la ley de Fourier: 


OT OT 
=-=k, A, —=-xk, dy da — 
A, Xx * Ox XxX ly a 
E E o. 
y Oy 

q, ==k, A, E =—k, dy — 


de donde: 


OX OX Xx 
O 
2 dy = -2- k,, dx dz ly a 2(x, TD lócayés 
Oy Oy O Oy Oy 
aa -2- k, dx dy jar 3 6 Jay dz 
Oz z Oz 


que sustituyendo de nuevo en la ecuación de la energía: 
2, ay ae pá k, ae dx dy dz [1 aye + q dx dy dz = A 
OX OX Oy Oy Oz Oz Ot 


Dividiendo por el volumen dx-dy-dz del elemento: 
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0 (%, —)+ aa O (% ra a (11) 
AN O O OO: Ot 


que constituye la ecuación diferencial general de la conducción del calor (o ecuación 
de difusión del calor) en coordenadas cartesianas, donde el significado de cada uno de 
los términos es el siguiente: 


(a) a Me 7 MUeéraz Potencia calorífica neta por unidad de volumen hacia el interior 


DEl” DE dx dy dz del elemento en la dirección £ (x, y Ó 2), [W/m']. 
q" Potencia calorífica generada internamente, [W/m'). 

ar e : 3 
ea Potencia calorífica almacenada por unidad de volumen, [W/m?]. 


En coordenadas cilíndricas la ecuación de conducción es: 


42 + dz 


E 
om 


A A 
4 N 
E 
+ 
E 


42” 
1 
k 
i 
F 
1 


4 [A dr +dr 


1 


S 
E 
Po 


y 


dz 


Volumen de control diferencial dr-r-d¿-dz para el análisis de conducción en 
coordenadas cilíndricas (r, y, Z). 


Eo (E): E E k, dE + e le lr E (1.2) 
r or or) r*00 0P ) Oz Oz Ot 
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Y en coordenadas esféricas: 


40 + de 


METTTZA  dór+dó 


3 


Volumen de control diferencial dr-r-sin 0-d4-d6 para el análisis de conducción en 
coordenadas esféricas (r, q, 0). 


E hor 1 > 1 E O k, OT + : 1 (o) [kasino 0) eq pc (1.3) 
r? or or) rísin*00pL *0p) r“sin0 00 00 ot 


Empleando el operador V, para materiales homogéneos (k,=k,=k,), la ecuación de 
conducción queda como: 


válida para cualquier sistemasde coordenadas elegido y donde: 


Coordenadas cartesianas: V=i—+j—+k 
Ox “0y Oz 


Coordenadas cilíndricas: V =i E a o + k E 
or "róop Oz 


Coordenadas esféricas: V =i S + ye k a E 
or “"róop  rsin0 04 


Algunos casos particulares interesantes de esta ecuación son: 


Material isótropo (k,=k,=k,=k) y conductividad constante (supuesto que no cambia 
con la temperatura): 
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Coordenadas cartesianas: 
O O ON 
AA O PO 


siendo q = EJ [m?/s] la difusividad térmica, otra propiedad de las sustancias. 
pc 


Coordenadas cilíndricas: 


rá ORO a On 
PORO on az O US a oe 


Coordenadas esféricas: 


Oo E A OE OIE OA 
a r E sin O o RS 
r?| or or) sin0 00 00) sin*0 06 k  ao0t 


Empleando el operador V?: 


mucho más sencilla, y válida para cualquier sistema de referencia, como era de 
esperar y donde: 


C. cartesianas: V? = 


Z 2 
C. cilíndricas: v?= nES l z ] a za + sl , 
rorkl or) r? op” Oz 
2 
C. esféricas: v?= 1130 G —)+ LE 30 (sino E |. e 50 
rólort. 0r) sing 00 00) sin*0 0p* 


Si además de tratarse de un material isótropo con conductividad constante, el 


ea a o, 
régimen es estacionario (2 = o): 
t 


Vr +7 =0 


Para material isótropo, conductividad constante o no constante (es indiferente), 
régimen estacionario, flujo unidireccional (por ejemplo, en la dirección x) y sin 


generación interna (q” = 0), empleando coordenadas cartesianas: 


¡04 y aE =0 = ERE at 
dxk dx dx 
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de donde: q, ES =-—k ar = cte (flujo de calor constante) 


m? dx 


Qs [w] =q,A= ae =cte (potencia calorífica total constante) 
x 


CONDICIONES DE CONTORNO Y CONDICIONES INICIALES 

Ecuación de conducción es: 

> De 1” orden en la variable temporal, exige una condición inicial: 
b=0 -T(EOJET, 


>» De 2% orden en las coordenadas espaciales, exige 2 condiciones de contorno en 
cada coordenada. Las condiciones de contorno habituales son: 


1. Condición de Dirichlet o de primer tipo: temperatura superficial constante (la 
superficie está en contacto con una sustancia en cambio de fase: 
fusión/solidificación o condensación/ebullición o con un fluido con el que 
intercambia calor por convección con un coeficiente de convección h muy 
elevado). 


Ts 


T(0,t)=Ts T«D 


2. Condición de Nuemann o de segundo tipo: flujo de calor constante. 


= Flujo de calor no nulo. Puede ser el caso de una superficie sobre la que 
incide una irradiación solar G; [W/m?]. q, correspondería con la fracción 


de energía radiante absorbida, esto es, q.,=0G,, siendo a la 
absortancia de la superficie. 


E T(<t) 


—E 
= Superficie aislada térmicamente (adiabática): 
T 
a OS 
Ola 
1x0 
—* 
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3. Convectiva o de tercer tipo: 


T(0) 3 
E KE Ide n.[r, —T(0,t)] ll ' STD) 
de 
4. Radiativa: : 
T(0,t) 
e peo T/ -T*(0,0)] To , 0%) 


También pueden darse condiciones de contorno mixtas en las que tengan lugar de 
forma simultánea, por ejemplo, intercambio radiativo con otra superficie, convección e 


irradiación solar (q, ): 


= so[r; —T*(0,£)]+ h.[T, —T(0,t)]+ q; 


x=0 


Sólamente cuando las condiciones de contorno son lineales (condiciones de Dirichlet, 
Neumann o convectiva con h¿ constante) pueden obtenerse soluciones analíticas a la 
ecuación de conducción. Esto no es posible con condiciones de contorno no lineales 
como es el caso de la condición radiativa o también de la condición convectiva con 


flujo laminar donde h, =h, (T, E yA En estos casos es necesario acudir a métodos 
numéricos de resolución de la ecuación diferencial. 
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2. CONDUCCIÓN ESTACIONARIA UNIDIMENSIONAL 


En este apartado se estudiará la conducción estacionaria unidimensional, 


? : : , : 0) 
caracterizada porque las variables objeto de estudio no varían en el tiempo (2 = o) 
t 


quedando definida la ley de temperatura en función de una única coordenada espacial. 
Tampoco se considerará, de momento, la generación interna (q =0). Los casos de 
estudio más significativos son los de paredes plana, cilíndrica y esférica. En todos 
ellos se supondrá materiales isótropos (la conductividad térmica no varía con la 
dirección considerada) y como aproximación que la conductividad del material es 
constante con la temperatura, lo cual es una buena aproximación a la realidad siempre 
que el rango de variación de temperatura dentro del material no sea muy amplio. 


Pared plana simple: sea una pared plana de espesor L, altura H, profundidad b 
(perpendicularmente al plano del dibujo) y temperaturas superficiales T, y T2. Para que 
la conducción a través de la pared sea  estricaamente unidimensional 
(perpendicularmente a la misma) es preciso que todo el canto perimetral esté aislado 
térmicamente. Si no dispone de tal aislamiento sólo se podrá asegurar la conducción 
unidimensional en las zonas alejadas del borde. 


Aislamiento térmico 


Ti Toa<T4 


TL 
L 
Xx 
Aplicando la ecuación diferencial de difusión del calor en cartesianas: 
0) oT 0) oT 0) oT € OT 
k + k, + k, +q =pce— 
OX Ox) Oy Oy ) Oz OZ Ot 
Que según lo establecido, se simplifica a: 
EA LO 
OXXk  OX 
y dado que sólo queda una variable independiente, la derivada parcial se puede 
expresar como total: 
El =S =0 
dx dx 


Y dado que la conductividad se ha establecido como constante: 


dr 


Integrando: de C 


1 
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T(x) =C,x+C, 


Las constantes de integración C, y C, se obtienen de la aplicación de las condiciones 
de contorno: 


T(x=0)=T1 T(x=L)=T2 
lo que da: 
Cy=-(T1-T2)/L Co=Tx 


Y por tanto: 


T(x) =1T, — (T, ple de (2.1) 


Ahora se puede utilizar la ley de Fourier para determinar la potencia calorífica: 


dT kA 1, =15 
=-—kA = T, —T, )- 
q a ñ (T, 2) 


Ka 
KA 


donde A=H-b y siendo = la resistencia de una pared plana. Expresando la ecuación 


(2) en forma gráfica: 


R=- 
Ta Ta 


De las expresiones de T(x) y q (ecuaciones 1 y 2), se deduce lo siguiente: 


"=  Latemperatura varía linealmente en la pared. 
= La potencia calorífica es constante a través de la pared, esto es, no depende 

de x. 

. T, -T sz Z 
=  Elflujo de calor | q = as e UE es también constante a través de la pared. 
A L/k 

=  Laresistencia térmica vale L/kA. 
"= Para una potencia q y un área A dados, el valor absoluto del gradiente de 
dr T-T 

E A A 

ye L kA 

decir, cuanto menor sea la conductividad térmica tanto mayor será la pendiente 
de la ley de temperatura. 


temperatura varía inversamente con la conductividad 


Otra forma de resolver la pared plana simple parte del hecho que q es constante, lo 
cual se podría haber deducido de la física del problema. Partiendo de la ley de Fourier: 


Separando variables e integrando entre la caras extremas de la pared: 
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Faax = —k Aj dT 


Resolviendo y ordenando: 


kA T-T 
q= L (T, ES e 


Para obtener la ley de temperatura, se puede integrar la ley de Fourier dejando 
indefinido el límite superior: 


Xx T(x) 
f, aax=-k Af] dT > T(9=1, Ex 


Llevando la expresión de q a T(x) resulta: 
X 
T(x) =1, — (T, -T, e 


Pared plana compuesta o múltiple: es la constituida por la superposición de varias 
paredes planas simples de espesores L; y conductividades k;. 


Tz Ta Ta Ta Tr+1 


Teniendo en cuenta que todas las paredes son atravesadas por la misma potencia 
calorífica: 


k,A 
Despejando las diferencias de temperatura y sumando: 


L 
T,-T, = E 
171) “LA 
za L, md 
TT = 9 O O: ia 
ZA 
L 
Led 


n 


Y en términos de resistencias térmicas: 
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L, 
Ca 


KA k A 
Tn+ ; 
ANA NN rd 


Son aplicables las observaciones hechas para la pared simple, pudiéndose comprobar 
cómo la pendiente de la distribución de temperatura en cada pared está en relación 
inversa al producto k-A. 


Pared cilíndrica simple: es la constituida por un tubo de longitud L cuyas caras 
interior (radio r,) y exterior (radio r,) están a temperaturas T, y T2, respectivamente. 
Para que la conducción a través de la pared sea estrictamente unidimensional (en la 
dirección del radio) es preciso que los extremos del conducto estén aislados 
térmicamente. Si no dispone de tal aislamiento sólo se podrá asegurar la conducción 
unidimensional en puntos no próximos a los extremos. 


Aplicando la ecuación de difusión del calor en coordenadas cilíndricas: 


10 or 10 o0Ti 0 oT s oT 
k,r AS k, + k, +q =pc=— 
r or or) r? 09 0p ) Oz Oz Ot 


Según lo establecido, se simplifica a: 


lr 20 
r or or 


y dado que sólo queda una variable independiente: 


Leto S) 20 
r dr dr 


Si la conductividad es constante, tal como se ha establecido: 
y 
=C 


Integrando: qe 1 
T(r) =C, Inr+C, 
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Las constantes de integración C, y C, se obtienen de la aplicación de las condiciones 
de contorno: 


T(r=r,)=T, T(r=r2)=T> 
lo que da: 
T, —T =T. 
E= 77 C, =1, += Inx, 
In + In + 
ñ ñ 
Y por tanto: 
Inr/r 
A a 2-3) 
Inr,/1, 


Ahora se puede utilizar la ley de Fourier para determinar la potencia calorífica: 


dT  27kL T, -T 
q =-kA ER (T, 1)e (2.4) 
dr mE nr) /r, 
n, 277 kL 
nr, /r, E dd 
donde A=211L y siendo a la resistencia térmica de una pared tubular. 
TT 


Expresando la ecuación (4) en forma gráfica: 


Ea a 
T, 277 kL 1, 


q 
De las expresiones de T(r) y q (ecuaciones 3 y 4), se deduce lo siguiente: 


"=  Lapotencia calorífica es constante a través de la pared. 


h LS as A 
= El flujo de calor |q = pe E disminuye para puntos tanto más hacia 
A inr,/r r 
el exterior del tubo. 
EOS nr, /r, 
"=  Laresistencia térmica vale —=—=. 
TT kL 


=  Latemperatura varía logarítmicamente con r: 


Ta>Ta Ty>T, 


Otra forma de resolver la pared cilíndrica parte del hecho que q es constante, lo cual 
se puede deducir de la percepción física del problema. Partiendo de la ley de Fourier: 
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q= ¿TE =-k(277 ma aa = cte 
dr dr 


Separando variables e integrando entre la caras extrema de la pared: 


Pa=aerLf, dr 


Resolviendo y ordenando: 
e T ES T, 
Inn /r, 
27 kL 


Para obtener la ley de temperatura, se puede integrar la ley de Fourier dejando 
indefinido el límite superior: 


In r/r, 


r dr T 
[[=-k2x 1[, ar A 


Pared cilíndrica compuesta o múltiple: es un tubo de varias capas, cada una 
caracterizada por su conductividad k; y sus radios interior r; y exterior r;,,, de tal forma 
que el exterior de una capa coincide con el interior de la siguiente. 


Operando análogamente a como se hizo con las paredes planas compuestas: 


q T, a Í, ni 


qe NR y y ran 


27 k¡L 


nm, /r, Inr,,, /r, 
27 k,L Ry 7 27 k.L 


T T Ta T Tn+ 
Ne VÁ y MA A ge To AA ro T m4 
q q 


Pared esférica simple: es la constituida por una esfera hueca de radio interior r;, y 
exterior r, y temperaturas Ty y T2, respectivamente. 


T 
q 
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ra 


Aplicando la ecuación de difusión del calor en coordenadas esféricas: 


271 qe + . za í 5 Le 0 [4 02 ed za 


r? Or or) r?sintgopl 2049) r?sino 00 00 Ot 
Que según lo establecido, se simplifica a: 
1 r? z =0 
r?* Or or 


y dado que sólo queda una variable independiente: 


1 d (kr) =0 
r? dr dr 


Y como la conductividad se ha establecido constante: 


eL 


Integrando: dh C, 


C 
T(r)=-=2+C, 
r 


Las constantes de integración C, y C, se obtienen de la aplicación de las condiciones 
de contorno: 


T(r=r,)=T, T(r=r2)=T> 
lo que da: 
_ 1-—T, a Li Lo 
Y, —Yr, A E 
Y por tanto: 


) Y/r, -Yr (2.5) 


A 


1 2 


Ahora se puede utilizar la ley de Fourier para determinar la potencia calorífica: 


dT 47 k T,-T 
a dr ME 2) Yr, —1/r, 8 
47 k 


Y/r, -1Y/r 
donde A=471? y siendo Mi téa la resistencia de una pared esférica. Expresando la 
TT 


ecvuación (6) en forma gráfica: 


p 1H V8 


7 Ark Te 


A A ro vu 


q 
De las expresiones de T(r) y q (ecuaciones 5 y 6), se deduce lo siguiente: 


=  Lapotencia calorífica es constante a través de la pared. 
= El flujo de calor |q = di ea disminuye para puntos tanto más 
A 1Yr, -Yr, r* 
hacia el exterior. 
Y/r, -Yr, 
dark 
=  Latemperatura varía inversamente con r/. 


=  Laresistencia térmica vale 


Ta>Ti Ti>T> 


Otra forma de resolver la pared cilíndrica parte del hecho que q es constante, lo cual 
se puede deducir de la percepción física del problema. Partiendo de la ley de Fourier: 


T T 
q= a E a —k(4xrr? ¡E = Cte 
dx dr 
Separando variables e integrando entre la caras extrema de la pared: 


f' E =-k4r [dr 


Resolviendo y ordenando: 


S Ls 
e 1 =U Ta 
477 k 


Para obtener la ley de temperatura, se puede integrar la ley de Fourier dejando 
indefinido el límite superior: 


Ear. T " Vr, -Yr 
[¡7=-kaz |, ar > E 


Pared esférica compuesta o múltiple: es la constituida por varias capas, cada una 
caracterizada por su conductividad k; y sus radios interior r, y exterior r;,,, de tal forma 
que el exterior de una capa coincide con el interior de la siguiente. Operando 
análogamente al resto de paredes compuestas: 
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e E am e E Ta 
A CEL 
E Pia 
47 k, 
121 LL 
R des Fi Pi = F; Pisa 
Ark, p 477 k; 
Ta3 T 
q q 


Existen una serie de geometrías en las que a pesar de no tener lugar conducción 
unidimensional (la ley de temperatura en su interior no se puede expresar en función 
de una sola coordenada), desde el punto de vista del cálculo de la potencia térmica 
que las atraviesa sí pueden tratarse como unidimensionales. Para ello es preciso que 
haya implicados sólo dos nodos térmicos (dos temperaturas límite diferentes) y exista 
resistencia térmica que caracterice la geometría. Es decir, que sea aplicable una 
expresión del tipo q=AT/R. Las configuraciones mostradas en la tabla siguiente, y en 
las condiciones expuestas, pueden tratarse como unidimensionales, presentándose 
sus resistencias térmicas. 
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Configuración 


Cilindros excéntricos 


Resistencia térmica 


2) co a e 
A > 
ZA = Es L>>E 
ns ES 27 kL 
P AS 
Prismas cuadrados concéntricos 0931 y — 0,0502 
R= z eo 
27 kL b 
0,7851n 2, 
R= b <1,4 
27 kL b 
L>>a 
Cilindro circular y prisma cuadrado 
concéntricos 
0,54 a 
ln === 
Ra==== M2 
27 kL 
Ma 
2h 
47 kr, 
La temperatura del medio en el infinito es T>. 
Cilindro enterrado cos) 
r 
R= 2 h<3n 
27 kL 
L>>r, 
2h 
2) 
ñ 
= h>3r, 
27 kL 


La temperatura del medio en el infinito es T>. 
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Configuración Resistencia térmica 
Viga rectangular enterrada 1 
R= ANO py 0078 
2,756 kL| In| 1+ — => 
a b 
L>>h,a,b 
La temperatura del medio en el infinito es T>. 
Arista con dos paredes adyacentes 

1 L 
R= =— W >= 

0,54kW 5 

W es la arista interna del cubo 
Esquina de tres paredes adyacentes 
T: 1 
a R=""— W>2 
¿LR 0,15k L 5 
ma Ls 
7, FIEL 
W 
Disco sobre un sólido semi-infinito 

SO 1 L 

e > 7 R=-—— W>-=— 

OZ sd ? 

La temperatura del medio en el infinito es T» 


Tabla 1. Resistencia térmica de paredes complejas (k representa la conductividad del material). 
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Cálculo aproximado de paredes compuestas heterogéneas. Sea una pared plana 
compuesta heterogénea como la de la figura: 


AAAAAAAAAS AAA 


Área A 


La Lac Lo 


Aunque el flujo de calor en algunos puntos del interior de la pared sea bidimensional 
q”=f(x,y), se puede obtener un aproximación de la potencia térmica total que atraviesa 
la pared bajo la suposición de flujo unidimensional, q”=f(x). Para ello caben dos 
modelos: 


e Las superficies perpendiculares a la dirección principal del flujo de calor se 
consideran isotermas, conocido también como el método de las rebanadas 


verticales. 


La 
q T; KA Ta 
» AAA 


e Las superficies paralelas a la dirección principal del flujo de calor se consideran 
adiabáticas, conocido como el método de las rebanadas horizontales. 


Con ambos circuitos térmicos equivalentes de paredes compuestas se obtienen 
diferentes resultados de potencia calorífica a través de la pared que se acentúan, 
discrepancia que se incrementa en la medida que las conductividades de las paredes 
en paralelo son tanto más dispares. Los valores obtenidos pueden considerarse como 
los extremos de un intervalo dentro del cual se encuentra la verdadera potencia 
térmica que a traviesa la pared. 
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3. CONDUCCIÓN ESTACIONARIA CON GENERACIÓN INTERNA DE CALOR 
UNIFORMEMENTE DISTRIBUIDA 


Se aborda en este capítulo el estudio de paredes con generación interna de calor q” 
[W/m*]. Tal es el caso cuando tiene lugar una reacción química exotérmica 
(desprendimiento de calor), una reacción nuclear, una corriente eléctrica o bien porque 
haya insertada una resistencia eléctrica uniformemente distribuida por el material. No 
hay que perder de vista que q” también puede resultar negativo cuando, por ejemplo, 
hay una reacción química endotérmica (absorción de calor). Se considerará en este 
capítulo que la generación de calor es uniforme en el material, lo cual responde bien a 
las situaciones habituales, y que la conductividad del material no cambia con la 
temperatura. 


Pared plana: considérese una pared plana de espesor L según la figura, cuyas 
temperaturas extremas T, y T, se suponen conocidas. 


q dv 
a E 
0) 
d, Arrdx = 4, me El dx a 4, as da, 
OX 
X dx 
AAA 
-1/2 : +L/2 
Figura 1.a. Pared plana Figura 1.b. Rebanada diferencial 


Partiendo de la ecuación diferencial de difusión del calor: 


0 OT 0 OT 0 OT m OT 
K, k k, Eq =pe— 
oxXL *0x) 0yL > 0y) 0” óz Ot 


Según las premisas establecidas, y como T=T(x): 


, ; A 
je =0 Ó k +q =0 


dx? 


Otra forma de obtener la ecuación diferencial hubiese sido considerar el equilibrio 
energético de una rebanada de espesor diferencial dx (figura 1.b): 


Centra ES U genera = d sale 


Y+ q Adx A da, 


Aplicando la ley de Fourier: 
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2 
dq, = E ae dx kad Pd 
dx dx dx 
Y sustituyendo en la ecuación anterior: 
2 
q Adx=-kA e a dx 
dx 
De donde: 
dEF 
k +q =0 
dx* Ñ 


Integrando dos veces resulta: 


” 


q 
EL)=S A +Cx+C, 


Aplicando las condiciones de contorno en los extremos de la pared. 
T(x=-L/2) =T, T=L2)=1 


se obtienen las constantes de integración: 


L, +T, 


ql ; X 
T(x) = L=LE)=+ 
0 MI 

AAA ]áZ- 


pared plana sin q'"" 


generación interna 


Que representa una parábola. Las diferentes situaciones que se pueden presentar 


son: 
T1=T max 
Ta 
2 Ta Ta 


Ley de temperatura con generación positiva (q>0). 


Ti Ta 
Ta 
To=T min 


Ley de temperatura con generación negativa (q<0). 


+ 
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Y la potencia calorífica q en un plano genérico x: 


1, -L 


7 q Ax + ñ 
generación interna KA 


pared plana sin q" 


Por tanto, las potencias caloríficas superficiales en las caras izquierda y derecha 
valdrán: 


L kA 


dT a 
del SAM TR 


=—k A 
U,1/2 dx 


x=-L/2 


Dado que la distribución de temperatura es parabólica, el punto de máxima/mínima 


Ed dT 
temperatura estará ubicado donde di = 0, para el cual: 
Xx 


T-D k 
XJer_, = 2 
dx L q 
De cuyo valor dependerá que el máximo/mínimo este ubicado dentro de la pared 
ó dl 
E ga o fuera Ec eA Sí El que la segunda deivadaS E Enrtlos sea 
L q 2 dx k 


negativa (q">0), (o positiva si q"<0), da la confirmación de que en x|er_, efectivamente 
dx 


hay un máximo (q>0) o un mínimo (q'"<0), de valor: 


T 


q q E 1 A(T, ES EA (T, —T,Fk ¿h +1, 

8k Lig” Eq 2 
Cuando q'>0, el flujo de calor a través de la cara fría (derecha en la figura 1.a) 
siempre tiene lugar hacia el exterior de la pared, mientras que en la cara caliente 
(izquierda figura 1.a) dicho flujo puede ser: 


= Hacia el interior —q" 2 +(T, ye IS para el cual T 
Z L La 2 
estará localizado exteriormente a la pared. 
L=B k 


= Hacia el exterior de +(T, 1) <0 C —=< z , en cuyo caso T' 


está dentro de la pared. 


=  Nulo —q. a +(T, py dla e a , estando situado T' justo en 
2 L Laa 


la cara caliente. 
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Si q <0, el flujo de calor en la cara caliente siempre es hacia el interior de la pared, 
mientras que en la cara fría dicho flujo puede ser hacia el exterior, hacia el interior o 
incluso nulo, dependiendo también de la situación de T', es decir, del valor del término 


dde frente a L/2. 
L q 


Como se ha podido comprobar, la localización de ds está ligada al sentido del flujo de 
calor en la pared caliente (q >0), o en la pared fría (q <0). En la tabla siguiente se 
recogen las diferentes posibilidades comentadas: 


TT k| L T-Tk|_L T-Tk|_L 
L ' 2 L q 2 L q 2 
T,=T 
qa >0 Qij2 (dijo > 0 5/2 
1/2 T, 
Ta Ta T 
q y <0 
q <0 dr 1/2 1/2 Quj2 =0 ZE > Sd 
Ta T¿=T* : 


Si las temperaturas de ambas caras coinciden (T,=T>2=T+), las ecuaciones se 
simplifican notablemente: 


8k e ] 

a09=q Ax 
A q a” >0 
A aa q <0 


En tal caso la distribución de temperatura es simétrica respecto al plano x=0 y, por 
tanto, la temperatura máxima/mínima está justo en el centro de la pared. Teniendo en 
cuenta las expresiones de Tmax!T min, T(X) Se puede expresar: 


q 2 " 
T(O)=T. -——x >0 
A) Tina >k q 
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do? " 
T(O)=T. —*—x <0 
(0) = Tin >k q 


Y las potencias caloríficas superficiales serán en ambas caras hacia el exterior (q>0), 
O hacia en interior (q <O): 


Bajo estas circunstancias (T,=T>=T,) es posible expresar la relación entre las 
temperaturas en el centro T(x=0) y la superficial T¿ con q mediante un modelo de 
resistencia térmicas. Particularizando q(x) en L/2: 


m L bd 2 1/2 
=q A— > AA 
412 =4 > q AT 
Y sustituyendo en las leyes de temperatura T max Y Tmin: 


ES a 
A 


">0 
max Ss 8k s AKA q 
is q1/2|L A 
E, AAC EPA q” <0 
8k AKA 
Y reordenando: 
mm E e == R = E 
q">0: d1/2 MS Tmax A AKA q 
AKA 1/2 
LET sen Am A 
q""<O: [Q,,2| == L == Ts Tinin 
4kA lay, 


Es decir, cuando T,=T>=T, es posible sustituir la mitad de la pared plana con 
generación interna por otra equivalente sin generación y cuya resistencia térmica es 


L 
a siendo L el espesor de la pared completa y no de la mitad. 
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Barra cilíndrica maciza 


Sea una barra cilíndrica maciza de radio R con generación interna uniforme q” y 
temperatura superficial Ts. 


Ey «e 


Urrdr q, El dr = = q, E dq, 


Figura 2.a. Barra cilíndrica Figura 2.b. Rebanada diferencial 


Partiendo de la ecuación diferencial de difusión del calor en coordenadas cilíndricas: 


10 oT 10 or. 0 oT 0 oT 
k,r ES k, + k, +q =pc= 
r or Or) r? 09 0p ) Oz Oz Ot 


Según las premisas establecidas, y como T=T(r): 


. 1d ( dr q 
lr «8 [re ).2 =0 


r dr dr 


Otra forma de obtener la ecuación diferencial hubiese sido considerar el equilibrio 
energético de una rebanada de espesor diferencial dr (figura 2.b): 


dentra de 4 genera = d sale 
q/+q 2xrLdr =q/+dq, (3.1) 


Aplicando la ley de Fourier: 


El rr EE jar = 22 kL A er (3.2) 
d dríi dr 


Llevando (2) a (1): 


q 27rL dr = -27kL Ar ir 
dri dr 


Simplificando y ordenando: 


dd re) + 
rdri dr 


En cualquier caso, integrando dos veces: 
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E 
T(r)= Til +C,Inr+c, 


= 0 (por simetría) y T(r =R)=T., 


a dT 
Aplicando las condiciones de contorno: añ 
r r=0 


” 


se obtienen las constantes de integración c, =0y c, =T, + ma La distribución de 


temperatura queda finalmente como: 


Y la potencia calorífica en un radio genérico: 


q(r) = tl =pLgr 
dr 


Que particularizada para el radio exterior de la barra: 


a(R)=q AR?L 


La temperatura máxima (q >0) o mínima (q <0) estará localizada necesariamente en el 
centro de la barra: 


” 


T = HL gs 
4k 


q ape 
Ak 


En la barra cilíndrica también es posible utilizar modelos de resistencias térmicas para 
representar la relación de la potencia calorífica q(R) que atraviesa la superficie de la 
barra con las temperaturas superficial T. y en el centro (Tmax O Tmin). Combinando las 
ecuaciones: 
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R= 
q>0: q(R) = ma la Tmax , 4TkL Ts 
47 kL q(r) 
p= 
o op T ed Toni 
— min 
q"<O: a(R)] a A s 3 
q) 
ArTkL 


Las ecuaciones de la barra cilíndrica maciza podrían haberse obtenido como un caso 
particular de otro más general constituido por una barra cilíndrica hueca de radios r, y 
r¿ y temperaturas superficiales T, y T2, como se muestra en la figura: 


Al aplicar las condiciones de contorno: 
T(r=r1)=T1 
T(r=r2)=T2 
a la ley general de temperatura: 


2 
T(r) e +c,Inr+c, 


resultan las constantes de integración c, y Cz: 


4k 
Qs > 
In + 
ñ 
q 
q rd r?) (T, 1) 
E =1, += Inr, 

ÑS In 
Fi 


La ley de temperatura queda: 


A e [O ao AO 


Por aplicación de la ley de Fourier, la potencia calorífica será: 
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2 2 
E 14 91 _ 27 kL(T, -T,) ó il 2 —r, má 


dr  In(r,/r) 21n(r, /r, 


que particularizada para r, y r,: 


E 27 kL(T, -T,) a TL e ] 
E pie caro 
_ 27rkL(T, —T,) "2 _1(r/r,Y | 
o e E 


Haciendo la diferencia entre las potencias caloríficas superficiales de la caras interna y 
externa: 


a(r,)-q(r,)=x Lq"(r? —rj) 


que, evidentemente, corresponde con la potencia generada internamente en la pared, 
dada por el producto de la generación interna q por el volumen de la misma 


mL (r? 12). 


Es posible establecer a priori el sentido de las potencias caloríficas en las caras 
interior y exterior del tubo por simple comparación de la situación del punto de 
pendiente nula en la distribución de temperatura (punto de máxima temperatura si 
q >0 o de mínima temperatura si q <0) con los radios r, y r,. Haciendo nula la derivada 
de la temperatura: 


resulta: 


Pr, 1 
a E 2 q 
r” =rjar_, === 
dr q In(r, /r,) 


siendo r* el radio donde dicha derivada es nula. La casuística que se presenta en este 
caso coincide con la que tiene lugar en una pared plana con generación interna. Es 
decir: 


= Con q >0 la potencia calorífica sale por la cara más fría, mientras que en la 
cara más caliente la potencia térmica puede ir hacia dentro de la pared, hacia 
fuera o ser nula. 


= Con q <0la potencia calorífica entra por la cara más caliente, mientras que en 


la cara más fría la potencia térmica puede ir hacia dentro de la pared, hacia 
fuera o ser nula. 
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